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摘 要 :复杂 工程 结构 中 往往 包含 多 源 不 确定 性 ,而 概率 -区 间 混 合 不 确定 性 理论 是 处 理 这 些 不 确 
定性 因素 的 有 效 手段 。 但 是 ,经 典 的 概率 -区 间 混 合 可 靠 性 计算 模型 由 双 层 优化 循环 吝 套 组 成 ,对 
于 复杂 工程 问题 存在 计算 成 本 过 高 和 鲁 棒 性 不 足 的 问题 。 为 了 改进 算法 性 能 ,提出 了 面向 概率 -区 
间 混 合 可 靠 性 分 析 的 自 适 应 共 力 控制 方法 ,该 方法 将 前 两 次 迭代 方向 与 当前 和 迭代 点 处 的 负 梯 度 方 
向 组 合 ,同时 加 入 控制 因子 对 组 合 后 的 方向 进行 修正 ,改善 其 迭代 点 的 下 降 方 向 ,提高 相应 收敛 速 
度 ,并 通过 可 行 下 降 方向 的 定义 证 明了 修正 后 方向 仍 为 优化 问题 的 下 降 方向 。 此 外 ,为 防止 迭代 过 
程 中 产生 数值 振荡 ,进一步 引入 了 Wolfe-Powell 准则 ,保证 控制 因子 自 适 应 le bes 
的 重 棒 性 。 为 了 获得 随机 参数 的 上 下 界 , 基 于 超 参 数 凸 模型 构造 了 相应 的 显 式 和 迭代 公式 。 计 算 

果 表 明 , 该 方法 在 解决 高 非 线性 问题 时 能 有 效 地 避免 迭代 过 程 中 i 
j 比 ,不仅 能 满足 精度 需求 ,而 且 可 以 高 效 、 稳 定 地 求解 可 靠 度 指标 。 
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Abstract: Complex engineering structures often contain multi-source uncertainties , while probability-inter- 
val hybrid uncertainty theory is an effective tool to deal with these uncertain parameters. However , the clas- 
sical probability-interval hybrid reliability calculation model is composed of double nested optimization 
loops , which has the problems of expensive calculation cost and insufficient robustness for complex engi- 
neering problems. In order to improve the performance, this study proposes an adaptive conjugate control 
method for probability-interval hybrid reliability analysis. This method combines the decrease directions of 
previous two iteration steps with the negative gradient direction at the current iteration step. A control fac- 


tor is also added to further modify the combined direction , aiming to improve the descent direction and 
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convergence speed. Besides ,it is proved that the revised direction is still the descent direction of the opti- 


mization problem according to the definition of the feasible descent direction. Moreover, in order to prevent 


the numerical oscillation in the iterative process, the Wolfe-Powell criterion is introduced to enhance the 


robustness of the algorithm by adaptive update of control factor. To obtain the upper and lower bounds of 


random parameters , the corresponding explicit iterative formulas are constructed based on the super para- 


metric convex model. The results show that the proposed method can effectively avoid the oscillation phe- 


nomenon in the iteration process when highly nonlinear problems are solved. Compared with the traditional 


method ,the proposed method not only meets the requirement for accuracy , but also solves the reliability in- 


dex efficiently and stably. 


Key words structural reliability; probability-interval hybrid uncertainty analysis; first order second mo- 


ment; conjugate direction method; adaptive conjugate control 


结构 工程 中 存在 的 大 量 多 源 不 确定 性 因素 会 对 
结构 安全 产生 严重 威胁 。 而 结构 可 靠 性 理论 可 
以 合理 地 处 理 这 些 多 源 不 确定 性 参数 ,保证 结构 的 
安全 性 。 一 般 来 说 ,不 确定 性 可 以 分 为 两 类 :随机 不 
确 起 性 和 认 知 不 确定 性 “ 。 随 机 不 确定 性 是 指 结 
构 趾 存在 的 固有 随机 性 , 它 不 能 通过 知识 或 者 数据 
的 增加 予以 消除 “1。 认 知 不 确定 性 源 于 知识 的 医 
过 或 者 简化 ,可 以 通过 获取 更 多 实验 数据 的 手段 来 
减 站 此 类 不 确定 性 91 。 

QJ 概率 方法 适用 于 处 理 随机 不 确定 性 ,其 核心 思 
柜 是 根据 大 量 的 实验 样本 数据 来 获得 结构 的 精确 可 
匈 麻 指标 ,经 典 的 计算 方法 包括 蒙特 卡 罗 模 拟 法 A 
方法 \ 代 理 模型 法 和 近似 可 靠 度 计算 方法 。 蒙 特 卡 
风 稿 拟 法 操作 简单 旦 精度 高 ,然而 其 高 额 计算 成 本 
严 量 阻碍 了 工程 应 用 。 为 了 减少 抽样 方法 的 计算 
量 6 学 者 们 相继 提出 了 多 种 改进 的 抽样 方法 ,如 重要 
TUERI UU .方向 抽样 法 "" 、 子 集 模拟 法 等 。 抵 
方法 的 求解 思路 是 通过 计算 结构 响应 的 高 阶 矩 来 获 
得 失效 概率 ,该 类 方法 求解 精度 较 高 ,但 受 样本 点 的 
位 置 影响 较 大 " 。 代 理 模型 法 使 用 了 显 式 数学 模 
型 对 真实 复杂 结构 进行 全 局 近似 , 大幅 提高 了 计算 
效率 ,但 该 类 方法 受 样本 点 的 影响 较 大 ,对 于 高 维 可 
靠 性 问题 的 适用 性 不 强 ' 5 。 近 似 可 靠 度 计算 方法 
中 代表 性 方法 是 HL-RF 和 迭代 算法 " ,该 算法 因 远 代 
公式 简单 .计算 成 本 低 而 被 广泛 应 用 于 可 靠 性 分 析 。 
但 对 于 高 非 线性 问题 , HL-RF 算法 容易 不 收敛。 对 
此 ,Yang ”提出 了 基于 混沌 控制 的 HL-RF 迭代 算 
法 ,以 实现 对 迭代 步 长 进行 稳定 控制 。 亢 战 等 ” 考 
虑 到 当前 迭代 点 和 它 在 极限 状态 曲面 内 法 线 的 夹 
角 , 并 通过 该 夹 角 来 修正 下 一 步 的 迭代 点 ,从 而 实现 
快速 稳定 收敛 。Gong 等 "提出 了 有 限 步 长 迭代 
法 ,通过 缩小 迭代 步 长 提高 了 可 靠 性 分 析 的 鲁 棒 性 。 


孟 增 等 分析 了 HL-RF 算法 解决 高 非 线性 问题 时 
发 散 的 原因 ,指出 了 迭代 点 振 萝 存在 方向 性 特征 ,并 
随后 提出 了 基于 修正 混沌 控制 的 HL-RF 算法 增强 
了 其 鲁 棒 性 。 

非 概率 方法 的 优势 在 于 对 实验 样本 数量 的 依赖 
METEO! ,因此 在 处 理 认 知 不 确定 问题 时 受到 学 者 
们 的 青睐 ,其 中 经 典 的 非 概率 模型 包括 区 间 模 型 和 
凸 模型 。 区 间 模 型 将 不 确定 变量 描述 为 区 间 形 式 ， 
通过 处 理 后 得 到 非 概率 可 靠 度 指 标 2 。 唐 嘉 旧 
等 .在 区 间 模 型 的 基础 上 建立 了 区 间 和 鲁 棒 优化 模 
型 ,有 效 地 处 理 了 带 有 不 确定 参数 的 非 线性 结构 鲁 
棒 性 优化 问题 。Wang 等 525 基于 区 间 模 型 对 超 音 速 
机 渐 结 构 进 行 了 非 概 率 可 靠 性 优化 设计 ,结果 表明 ， 
安全 系数 法 在 描述 设计 过 程 中 的 不 确定 性 方面 存在 
缺陷 , 且 区 间 模 型 可 以 精确 地 表征 少量 实验 样本 的 
不 确定 性 参数 ,更 接近 实际 工程 问题 。Gurav 457 
发 现 区 间 模 型 给 出 的 最 危险 情况 是 取 每 个 变量 的 极 
值 。 江 涛 等 ' 吕 采用 反射 运算 技术 求解 了 区 间 模 型 
的 非 概率 可 靠 度 指标 。 相 比 于 区 间 模 型 , 凸 模型 则 
侧重 于 描述 随机 参数 间 的 相关 性 '” ,因此 近年 来 受 
到 了 广泛 关注 。Canzerli 等 中 基于 凸 模型 提出 了 一 
种 不 确定 荷载 下 的 结构 优化 设计 方法 。 亢 战 等 1 
认为 基于 凸 模型 的 最 优 拓扑 构 型 设计 比 确定 性 的 拓 
扑 构 型 设计 更 为 安全 有 效 。 由 于 可 供 选 择 的 经 典 非 
概率 模型 过 少 , Meng 等 ”构造 出 更 广义 的 超 参数 
凸 模型 ,实现 了 有 限 实验 样本 的 高 置信 表征 。 

实际 工程 中 往往 同时 包含 随机 和 认 知 不 确定 
性 ,因此 学 者 们 把 概率 模型 和 非 概 率 模型 结合 后 提 
出 了 混合 不 确定 性 模型 。Elishakoff 等 1 深入 研究 
了 概率 模型 和 凸 模型 的 混合 不 确定 性 问题 。Qiu 
等 六 将 概率 方法 与 区 间 分 析 结 合 , 并 通过 两 个 析 相 
系统 证 明了 该 方法 的 可 行 性 。Kang 等 将 多 源 不 
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息 来 获得 精确 的 随机 参数 ,相应 结构 的 
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融合 后 提出 了 混合 可 靠 性 建 模 方法 。 刘 海 波 等 RIR EUT 

D ro 不 确定 性 参数 的 系统 可 靠 = PIG) x 01 

性 分 析 方法 ,并 指出 该 方法 只 需要 少量 的 实验 数据 - [caf GO dr 

即 可 确保 系统 安全 ,适合 于 复杂 工程 结构 。Jiang 

等 "进一步 指出 ,在 实验 数据 不 足 的 情况 下 ,一 些 = ed (1) 

重要 的 分 布 参 数 ( 如 均值 .标准 差 等 ) 无 法 直接 确定 ”其中 :xz = [ae e MARAA n 维 随机 参数 ; 

而 是 以 区 间 的 形式 存在 。 此 时 ,概率 -区 间 混 合 不 确 AOO Fel ) 分 别 表 示 功 能 函数 G(x) 的 概率 密 
定性 模型 是 处 理 这 类 问题 的 有 效 手段 。 BE PRHBURI STE A f PC, i C(x) >0 时 表示 结构 安 
尽管 概率 -区 间 混 合 不 确定 性 模型 比 概率 和 非 ” 全。 实际 工程 中 往往 难以 直接 求解 上 述 高 维 积 4 

概率 模型 更 适合 于 处 理 复杂 的 工程 问题 ,但 是 a 因此 通常 采用 一 次 可 靠 度 方法 来 计算 失效 概率 和 可 


性 能 仍然 不 足 。 概 率 - 区 间 混 合 不 确定 性 模型 通 
et ey 
对 于 高 非 线 性 函数 ,概率 -区 间 混 合 不 确定 性 模型 存 
在 延 代 振荡 和 收敛 性 问题 。 因 此 ,起 待 发 展 一 种 高 

性 能 的 计算 方法 ,在 增强 混合 模型 收敛 性 的 同时 降 
ORE. 

加 为 了 实现 上 述 目标 ,本 研究 提出 了 基于 自 适应 
共 殉 控制 的 概率 -区 间 混合 可 靠 性 计算 方法 。 首 先 ， 
握 尖 了 自 适应 共 罗 控制 方法 ,其 核心 思想 是 将 前 两 
次 枉 迭 代 方 向 与 当前 先 代 点 处 的 负 梯 度 方向 组 合 形 
Wife E (43588 75 168] ; 其 次 ,把 该 方法 拓展 至 概率 -区 
呈 渴 合 可 靠 性 计算 中 来 控制 其 交代 的 方向 ,保证 算 
法 铺 棒 性 的 同时 提高 其 计算 效率 ;最 后 ,为 了 获得 随 
机 委 数 的 上 下 界 ,基于 超 参数 上 模型 构造 了 相应 显 
sepes 


1 三 概率 模型 区间 模 型 与 概率 区间 
'O 混合 不 确定 性 模型 


1.1 概率 模型 


概率 模型 适用 于 处 理 随 机 不 确定 性 ,其 需要 大 


O [A 
(a) 收敛 


图 1 


徘 度 指标 6。 在 进行 可 靠 性 分 析 之 前 ,需要 通过 
Rosenblatt 变换 将 随机 变量 x 从 原始 空间 转换 到 标 
准 正 态 空间 , 则 功能 函数 从 G(x) 转 到 g (uw) ,其 中 
gCu)d&m G(x) 在 标准 正 态 空 间 下 的 功能 函数 。 因 
此 ,B 的 几何 意义 就 是 在 标准 正 态 空间 中 坐标 原点 
到 极限 状态 曲面 的 最 小 距离 。 其 数值 计算 可 归结 > 
以 下 优化 问题 。 


ming- |u | 
d (2) 


s.t. g(u) 20 
上 述 优化 问题 常用 HL-RF 算法 求解 ,其 迭代 公 
式 如 下 。 


k+l 


kl _B Vg (uy a) 


HLRF = (3) 
id | Vg (Uti ai) | 
Ba g( Uire) -Vg Dru Winds (4) 
|l Vg (Uir) | 


其 中 ,Vg( Uher) XZR g (whar) 的 梯度 向 量 。HL-RF 
算法 迭代 公式 简单 .计算 成 本 低 , 但 对 于 强 非 线 性 问 
Bi TEENS, WWE 1 所 示 , HL-R 算法 在 二 维 空 
间 中 会 出 现 不 收敛 问题 。 因 此 ,学 者 们 提出 了 不 
同 的 改进 方法 “| ,但 这 些 改 进 算法 都 只 考虑 了 随机 
不 确定 性 ,没有 考虑 到 混合 不 确定 性 问题 。 


(b) 不 收敛 


两 种 典型 的 迭代 历史 


Fig.1 Two typical iteration histories 


投稿 网 站 :http://ejam. xjtu. edu. en ” 微 信 公众 号 :应 


] 力 学 学 报 


& 6 期 


1.2 区 间 模 型 


非 概率 方法 因为 对 实验 样本 不 敏感 而 常 被 用 于 
处 理 认 知 不 确定 性 参数 。 区 间 模 型 是 常见 的 非 概 率 
模型 ,其 计算 公式 可 以 表示 为 


min G(x(1 ,0)) 
| (5) 


s.t. xX srar" 
为 方便 计算 ,可 以 将 式 (5 ) 的 区 间 参 数 用 q; = 
(x,-2,,)/2,(121,2,3, ^, n) VEIT EHE, x FE xo 
: HS 的 中 点 和 半径 , 则 式 (5) 转 


(6) 
LABES! 


HE :4 是 归 一 化 的 不 确定 变量 ;g(9) 表 示 变 量 转 换 
局 的 功能 函数 。 对 于 4 空间 中 的 二 维 问 题 ,其 形状 
O 个 边 长 为 2 的 正方 形 , 如 图 2 所 示 。 


E 


s.t. max{ |q l, lgl, lg l; 


&q)70 


&(q)-0 


——— ERE EE 


rz 图 2 二 维 区 间 模型 形状 图 


l € Rig. 2 The two-dimensional interval model shape graph 
1.3 概率 -区 间 混 合 不 确定 性 模型 


概率 -区 间 混 合 不 确定 性 模型 综合 了 概率 和 区 
间 模 型 的 优势 而 被 用 于 处 理 复 杂工 程 中 的 混合 不 确 
定性 参数 。 该 模型 中 涉及 到 内 外 两 层 座 套 的 优化 循 
环 ” :外 层 循环 求解 可 靠 度 指标 的 上 下 界 ,内 层 循 
环 根据 区 间 参 数 寻 找到 功能 函数 的 上 下 界 。 其 中 ， 
可 靠 度 指标 下 界 B 的 求解 模型 为 


B - min ||u | 
7 
et =0 (7) 
G(x(u ,m)) =min G(x(u ,0)) 
| L U (8) 
s.t x <x<x 


$6] TY, S, MARRERO TAE E OLEUM 


可 靠 度 指标 上 界 B 的 求解 模型 为 


B - min || u | 
: (9) 

s.t. g(u) =0 
G(x (ji ,0)) =max G(x(p ,0)) (10) 


U 
xzx 


E Jr H jp f h AS AE E- X E] 


2 基于 自 适 
混合 可 靠 性 计算 方法 


2.1 AEEA A 


4189875 T5] A A ATARA rb FI fp E 75 1] 
5 EPRI T8) 28 JE JUS] P EET In], SCR 
快速 收敛 的 效果 "|。 但 是 对 于 混合 可 靠 性 计算 问 
题 , 共 红 方向 不 一 定 为 最 佳 下 降 方向 , 且 由 于 可 靠 度 
KORIRA E B 的 径 向 ,对 于 整体 步 长 控制 会 带 来 额 
外 计算 成 本 。 因 此 本 研究 提出 自 适 应 共 力 控制 方 
法 。 如 图 3 所 示 ,CD KRYEN HY , AC 表示 自 适 
mdtgefedrik. Bas wtf 7715 Hj BU AAKE 
代 方 向 对 当前 迭代 方向 进行 修正 ,并 加 入 控制 因子 
对 其 进行 控制 ,提高 算法 的 鲁 棒 性 和 计算 效率 。 此 
外 ,由 于 振荡 方向 沿 B 的 径 向 ,可 以 通过 (Bd *)/ 
| 人 ”| 保证 其 在 径 向 上 B 的 长 度 。 因 此 , 自 适应 
共 斩 控 制 方法 的 公式 可 表示 为 


- ; d'* 
usc =B | d | , 
k+l Vg" Qd d 
m a TER 
jas 
qa» 
g ECCE wo (11) 
EZ 


Hep ue Adsl DSL T og SEHE 8 75 T E E 3 FE FB ; 
为 防止 上 一 步 迭 代 方 向 对 当前 迭代 方向 修正 过 当 ， 
& 的 上 限 设置 为 0. 5; g! H guio); Vg 为 
Veg( uac)o 

当 功 能 函数 振荡 时 , 即 e^ -gf >0。 为 找到 最 
优 下 降 方 向 ,可 以 引入 Wolfe-Powell 准则 自 适应 控 
dl 大 小 。 具 体 表示 如 下 。 
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fne) Ol) uUctc(g), 
0. 5x f( u^, +a" d) f(u^.) *o'p Vf (uw. )' d 
K =I 2KVf (us, to d') d'2o Vf (us) d. 
K' else 
(12) 
Hop uut BERG o^ 是 迭代 步 长 , 即 为 Bi;c = 
10? ;0«p«li;p«o«l, 


CN 图 3 mgasw3tgefsml EXER SR 
- 3. The iterative history of "ud conjugate control method 
N 引 理 : 如 果 满 足 c>y,x* »o,H 60K o 
lupus 04) « (13), M 2 39 f CO BET 
KE, EIV (ust  )'d**' <0, 
证 明 :根据 式 (12) , 罚 函 数 对 uic 的 灵敏 度 为 


| VfC uic ) 22 us; +2cg Vg" (13) 
= Vg" | d' , odo 
= — +w E " -0 (14) 
Evye aeg? qae 
i (uic) *ug (Vg) =0 (15) 


RHE F EE 77 I8] HJ x X, 24 (VfCu c ))'d**' «0 
HF, d' B fu FE RFH. 
(VfCus d" = 


| Vg : d' à d^ 
2 (u! "( zi ppp i 
Gn) ve “len fal 
| l Vg : d! 7 d^ 
2cg^ (V Y - «e +ô = J 
SONIS | vel lel “~ aw 
(16) 


考虑 到 最 可 能 失效 点 处 e" =0。 因 此 , 当 g*>0 
Ih] e" 必须 要 下 降 ,此 时 (Ve*)"d* <0; g^ «O 时， 
& 必须 要 上 升 ,此 时 (Ve*) d' >0。 所 以 
g ( Vg^)! d' «0 
由 于 d' Jg E kA ISI B. "0 , 则 


(17) 


应 用 力学 
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d! 

sign| g" (Vg")" )= 

e ja 

. . d ' dt~! 

sign( g" SZAH | d | +ô | d | J 

AIX kf >w, PF 

, ; Vo" 4( d , d! 
‘(Ve ' -8 «e| — 46 IE 
s (V8) | vtta qae 

; V pa , od 

“(Ve i UE e| ~ +ô z J =0 
€ (e) | -Tyr eg? jal 

(18) 


此 时 ,车 c>y,x* >w, 且 6 一 0。 将 式 (18) 带 入 
到 式 (16) 中 ,得 
( VfCcust ))'d* = 
Vg“ " d y dt~! 
ive Tr Te 
Vg" " d! E a -1 
ves tte * ja71)]* 
Vg k d EL dU 
ivey tert ia71)])* 
k kT Ver k d 5 d 
ye! (Ve) "| ed Uer * je71)] = 
(19) 
综 上 所 述 ,当心 >w,c> 且 8 一 0 时， 
( VfL uz ))'d'*' <0, 说 明 dr*! 为 可 靠 性 指标 的 下 
降 方向 ,同时 也 证 明了 自 适应 共 思 控 制 方 法 的 可 
行 性 。 
对 于 概率 -区 间 混 合 不 确定 性 模型 ,其 外 层 优化 
可 以 通过 自 适应 共 轿 控制 方法 计算 ,对 内 层 优化 ,以 
求解 B 为 例 ,考虑 到 max| 1g11,1g,1,1g31,…,1g,1) = 
lim [|a |, 1, Lg 的 求解 问题 可 以 转换 为 超 参数 
凸 模型 的 求解 问题 , 即 


Ce 


2cg ( ve 


2 Geo - 


min £(q) 
|. t. lim|g|,«l Dd 
其 中 :4 是 非 概率 空间 中 归 一 化 的 不 确定 性 变量 ;p 


需要 设置 为 足够 大 的 正 数 以 保证 计算 精度 ,本 研究 
p 741000, HX (20) Karush-Kuhn-Tucker 最 优 条 
件 , 可 以 推导 出 以 下 迭代 公式 二 。 


dg(g) = 
= 9g(q) ðq; E NN 
E sien| ðq; ) e | Lemka] 


(21) 
其 中 表示 随机 参数 的 数量 。 同 理 , 模 型 内 层 优化 
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迭代 式 为 
ag(q) 1 
2 sion| 22D ðq; gs] 
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2.2 灵敏 度 分 析 


在 可 靠 性 计算 中 灵敏 度 分 析 必 不 可 少 , 使 用 有 

限 差分 法 可 以 直接 获得 3g(q")/0q;。 但 对 于 涉及 有 

限 元 分 析 的 复杂 工程 结构 ,有 限 差分 法 的 计算 成 本 

比较 大 ,此 时 可 以 用 摄 动 法 进行 有 限 元 灵敏 度 分 
析 '”]。 对 于 静 力 问题 ,结构 平衡 方程 可 表述 为 

K(q)U(q) -f(q) (23) 

"HE q = qo + Aq 由 一 个 不 确定 变量 gq。 和 一 个 

摄 动 Aq 组 成 。K e R" ”表示 整体 刚度 矩 阵 ， 

UeR'" 表示 节点 位 移 向 量 ,fe R 表示 外 部 和 荷载 


pr K AI U 对 go 的 一 阶 泰勒 展开 式 分 别 为 
K(q) = Kq) + >» XT (24) 
-一 oU 
N UCg) = Uq) + > 9g, Mi (25) 


eS EROR) ~ CA 分 别 考虑 零 阶 
无 穷 小 量 和 一 阶 无 穷 小 量 ,可 得 到 9U/9g;。 根 据 链 
式 法 则 ， ag( q^ )/aq: 可 以 通过 以 下 公式 计算 。 


9g(q) - g") oU 
ec ðq; ðU ag (26) 
Cji 
= üg(u') _əg(u*) , ðU 
Q du 0U gu Veg 
此 外 ,也 能 采用 摄 动 法 推导 获得 应 力 约 束 的 灵 
ARE, 


2.3 计算 流程 


本 研究 所 提 方 法 的 计算 流程 图 如 图 4 所 示 , 计 
算 步 又 概述 如 下 。 

1) 定 义 功能 函数 g(x)。 设 置 最 大 迭代 次 数 为 
200,k 20,0 2600,72 210 5, ^ =1。 

2) 将 随机 变量 归 一 化 处 理 得 到 9。 

3) 通 过 式 (21) ~ (22) 计算 归 一 化 随机 变量 q 
的 上 下 界 。 

4) 通 过 Rosenblatt 变换 将 q 的 上 下 界 转化 到 标 
准 正 态 空间 。 

5 ) 根 据 式 (11) 计 算 可 靠 度 指 标 上 下 界 。 
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6) 若 迭代 过 程 中 产生 振荡 , 则 使 用 Wolfe-Powell 
准则 自 适应 更 新 e 的 大 小 。 

7) 计算 前 后 两 次 迭代 点 之 间 的 误差 | x^ -x 7 
| x] 入 se。 如 果 满 足 误差 要 求 则 停止 迭代 ,如 果 不 
满足 则 转 至 步骤 (2 ) 重复 计算 直至 满足 收 剑 要 求 。 


定义 功能 函数 ， 设 置 初始 值 


随机 变量 


A a 
量 区 间 g 的 上 下 


通过 Rosenblatt 变 换 将 


q 
的 上 下 界 转 化 到 标准 正 态 空 i 


根据 式 (11) 计 算 可 靠 度 指标 上 下 界 


若 达 代 过 程 中 产生 振荡 ， 则 使 用 
Wolfe-Powell 准 则 自 适应 控制 x 的 大 小 


图 4 基于 自 适 应 共 轿 控制 的 概率 -区 间 
混合 可 靠 性 计算 方法 
Fig.4 The probability-interval hybrid uncertainty reliability 


calculation method based on adaptive conjugate control 


3 算 H 


在 本 节 中 ,通过 一 个 数学 算 例 、 一 个 矩形 悬臂 梁 
算 例 和 一 个 工业 机 右 人 应 用 来 验证 本 研究 所 提 方 法 
的 效率 和 精度 ,并 分 别 与 序列 二 次 规划 ( sequential 
quadratic programming,SQP) 算 法 、HL-RF 算法 混沌 
控制 (chaos control, CC) Zr E RIED HEUS XE VETE HRS 
方法 (chaotic conjugate stability transformation meth- 
od, CCSTM) 进行 对 比 。 所 有 算法 收敛 精度 为 | x - 
x^ Z xt 1 和 10…; 最 大 迭代 次 数 200;A =0. 1; 
K =1;N 表 示 性 能 函数 调用 次 数 。 对 于 自 适应 共 
斩 控 制 方法 ,矩形 悬臂 梁 算 例 的 灵敏 度 计 算 方 式 为 

最 动 法 ,其 余 均 采用 有 限 差 分 法 。 
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3.1 二 维 非 线 性 功能 函数 


第 一 个 算 例 是 一 个 二 维 非 线性 功能 函数 ,随机 
变量 Xi ~ N (ua ,0,4) 和 X5 ~N(jw,0w) 是 服从 正 态 
分 布 的 两 个 随机 参数 。 js = [mi -0. lo; p+ 
0.1o;],0,; 2[o;-0.1o;,0; *0.1o;],i 2 1,2 是 随 
机 参数 的 波动 区 间 , 其 中 (mw uus) = (10,9.9), 
(01,02) =(5,5)。 计算 B 下 界 和 上 界 的 初始 迭代 
点 分 别 设置 为 1 -0.10; 和 +0. 10;。 功 能 函数 定 
义 为 

G(x) xx) *x, x, * x, —18 (28) 

d 1 中 分 别 列 出 了 5 种 不 同方 法 的 计算 结果 。 
由 表 可 知 ,HL-RF £k fonti Hi f PE DR 77 E 
在 迭代 过 程 中 会 出 现 数值 不 稳定 现象 ,导致 迭代 无 
法 收敛 。 混 沌 控制 方法 达到 最 大 迭代 次 数 仍 未 收敛 
到 最 可 能 失效 点 。SQP 算法 能 稳定 收敛 到 最 可 能 
效 点 ,计算 结果 也 与 本 研究 所 提 方 法 相同 ,但 其 函数 
调用 次 数 较 多 。 尽 管 自 适应 共 红 控制 方法 在 求解 过 
各 中 同样 出 现 过 数值 不 稳定 现象 ,但 利用 自 适应 共 
固 弄 降 方 向 克服 了 收敛 性 困难 的 问题 。 结 果 表 明 
Walfe-Powell 准则 能 自 适应 控制 « 的 大 小 从 而 找到 
更 沽 适 的 迭代 方向 ,在 保证 算法 稳健 性 的 同时 提高 
IARR, R2 中 展示 出 在 不 同 初始 迭代 点 下 本 
研究 方法 的 计算 结果 。 从 表 中 可 以 看 出 自 适应 共 罗 
控制 方法 对 初始 点 的 位 置 不 敏感 。 

A1 算 例 1 各 方法 计算 结果 


Tab.1 The calculation results of each method 
= of the first example 
FEE B Na B Na 
SQP 2.218 108 2.397 96 
HL-RF - = = - 
CC - - - - 
CCSTM - = 三 z 
Proposed 2.218 80 2.397 68 


表 2 算 例 1 不 同 初始 点 计算 结果 
Tab.2 The calculation results of different 


initial points of the first example 


初始 点 B Noa B Na 
[5,4.9] 2.218 75 2.396 82 
[10,9.9] 2.218 80 2.397 67 
[15,14.9] 2.218 86 2.397 66 
[20,19.9] 2.218 86 2.397 66 
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本 算 例 考虑 一 个 三 维 悬 臂 梁 ,并 通过 有 限 元 分 
析 获 取 结 构 响 应 。 如 图 5 所 示 , 梁 的 长 为 640 mm, 
宽 为 200 mm ,厚度 为 50 mm , 左 端 固定 , 右 端 受 不 确 
定 荷载 。 甚 辟 梁 的 最 大 允许 位 移 不 超过 3. 8 mm, 
表 3 中 列 出 了 梁 各 随机 参数 的 大 小 和 分 布 ,随机 参 
数 波动 区 间 为 ws。 = [us -0. 10; u; *0.10;], 0, = 
[c; -0.10,,0;+0.10;],i=1,2,3。(w; -0.1c;) 和 
(m: +0. 10;) 分 别 为 计算 B 下 界 和 上 界 的 初始 点 。 
功能 函数 定义 为 ”| 

G 23.8 - A(E,v) (29) 
其 中 :5 是 杨 氏 模 量 ;v 是 泊 松 比 。 

计算 结果 如 表 4 所 示 。 可 以 看 出 所 有 方法 都 能 
稳定 收敛 , 且 HL-RF 算法 混沌 控制 方法 和 混沌 共 
力 稳 定性 转换 方法 算出 的 B 区 间 与 本 研究 方法 计算 
结果 相同 。 但 本 研究 方法 的 计算 效率 约 为 HL-RF、 
混沌 控制 和 混沌 共 斩 稳 定性 算法 的 1.5 倍 16 倍 和 
14 倍 。 此 外 , 相 较 于 SQP 算法 ,本 研究 方法 的 效率 
提高 了 11 倍 以 上 , 且 可 以 计算 得 到 最 佳 的 B 区 间 。 
尽管 本 研究 方法 在 计算 过 程 中 会 出 现 数值 振荡 现 
象 ,但 通过 自 适应 更 新 e 的 数值 解决 了 该 问题 。 
此 , 自 适应 共 斩 控 制 方法 对 于 本 算 例 是 非常 有 效 和 
精确 的 。 


640 mm 


p mm 


网 mm 
图 5 EESE TRA ASIE 


Fig.5 The structural diagram of rectangular 


cantilever beam 
表 3 PERBARISAN fs 
Tab.3 Random parameter distribution of 


rectangular cantilever beam 


参数 i Hi gi 分 布 
E/GPa 1 210 5 正 态 
v 2 0.3 0.01 正 态 
F/kN 3 200 20 正 态 
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重 增 ， 等 : 面向 概率 -区 间 混 
AKA 算 例 2 各 方法 计算 结果 


Tab.4 The calculation results of each method 


of the second example 


算法 B Nean B N oan 
SQP 2.881 288 3.393 340 
HL-RF 2.765 36 3.529 44 
CC 2.765 452 3:529 460 
CCSTM 2.765 408 3.528 408 
Proposed 2.765 26 3.529 28 


3.3 六 自由 度 工业 串联 机 器 人 


本 算 例 是 一 个 六 自由 度 工 业 串联 机 器 人 , 它 通 
一 系列 连 杆 相互 连接 ,如 图 6 ~7 所 示 。 机 器 人 未 

行 器 的 位 置 精 度 是 机 器 人 重要 的 性 能 指标 之 
-© 其 定义 为 实际 位 置 (mi ,ms m.) 与 末端 执行 器 所 
TAEC ,n,n) 的 偏差 。 表 5 中 列 出 了 机 器 人 各 
随机 参数 的 大 小 和 分 布 。 每 个 连 杆 的 长 度 与 前 三 个 
状语 角度 被 定义 为 服从 极 值 工分 布 的 随机 参数 。 由 
开间 三 个 关节 角度 对 位 置 可 靠 性 的 影响 较 小 ,在 此 
被 守 义 为 确定 性 变量 。 na = [jp; -0. 10;, p; + 
0:160,;],0,;7[o; -0.1o;,0; +0. 1o9;],i 2 1,2,3, 
DO 为 随机 参数 的 任 一 波动 区 间 。 计 算 B 下 界 和 
上 上 注 的 初始 点 分 别 设置 为 (w -0 lz) 和 (Am + 
0-le;). uf RC 7 

E nc + (m, =n,)? + (miemjg -0. 4° 


- (30) 
t " 
Orza HXH 轴 J4 fih J5 hh J6 
连 杆 4 HERE 6 
发 动机 5 
X — a. 
发 动机 4 mm T M 
DRE 3 Bs 
腕 关节 
Ly: 
XH 
动机 3 
RA 发 动机 
Bets 
发 动 
到 6 六 自由 度 串 联 工 业 机 器 人 结构 示意 图 


Fig.6 The structure diagram of 6-DOF series industrial robot 
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Put 
L 


Fig.7 The schematic diagram of the coordinate 
system of the ith connecting rod 
表 6 中 展示 了 5 种 不 同方 法 的 计算 结果 。 如 表 

所 示 ,SQP 算法 不 能 完整 地 计算 得 到 B 区 间 , 这 表明 
对 于 复杂 的 功能 函数 ,SQP 算法 的 鲁 棒 性 相对 较 弱 。 
HL-RF 算法 混沌 控制 方法 混沌 共 斩 稳定 性 转换 方 
法 和 基于 自 适应 共 红 控制 方法 具有 更 好 的 鲁 棱 性 
与 HL-RF 算法 和 混沌 控制 方法 比 , 自 适 应 deed 
方法 的 计算 效率 明显 提高 。 混 沌 共 轿 稳定 性 转换 方 
法 的 计算 结果 与 其 他 方法 相 比 误差 较 大 , 且 该 方法 
的 计算 效率 过 低 。 因 此 , 自 适 应 共 斩 控 制 方法 具有 
非常 明显 的 效率 和 和 鲁 棒 性 优势 ,也 为 处 理 混 合 不 确 
定性 问题 提供 了 有 效 的 计算 工具 。 

表 5 六 自由 度 串联 工业 机 器 人 随机 参数 分 布 

Tab.5 Random parameter distribution of 6-DOF 


series industrial robot 


参数 i Hi gi VES 
L,/mm 1 585 0.050 Ti fé T 
L,/mm 2 181 0.055 IKIE I 
L,/mm 3 665 0.015 极 值 
L,/mm 4 724.5 0.014 极 值 I 
Ls/mm 5 0 0.050 极 值 
Le/ mm 6 90 0.010 极 值 
01/(°) 7 -14.5 5.7x107^ 极 值 I 
6,/(?) 8 34.6 1.5x102 — AI 
6,/(?) 9 -10.2 2.9x102 极 值 I 


表 6 算 例 3 各 方法 计算 结果 
Tab.6 The calculation results of each 


method of the third example 


算法 B Nan B Nan 
SQP - - 3.037 480 
HL-RF 1.531 130 3.074 160 
CC 1.531  1.220x10 4.074 1.270 x10? 
CCSTM 1.897  L011x10? 3.826 1.011 x 10? 
Proposed 1.532 62 4.075 82 
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本 研究 针对 实际 工程 中 普遍 存在 的 混合 不 确定 
性 问题 展开 研究 。 首 先 ,在 共 轿 方向 法 的 启发 下 , 提 
出 了 自 适应 共 轿 控制 方法 ,并 将 其 应 用 到 概率 -区 间 
混合 可 靠 性 分 析 中 。 自 适应 共 罗 控制 方法 考虑 了 大 
代 点 振荡 时 的 方向 性 规律 ,将 前 两 次 迭代 方向 与 当 
前 迭代 点 处 的 负 梯 度 方向 组 合 ,并 加 入 了 控制 因子 
将 组 合 后 的 方向 进行 修正 ,从 而 大 幅 提高 了 计算 效 
率 。 同 时 ,根据 可 行 下 降 方向 的 定义 ,证 明了 修正 后 
的 方向 为 可 靠 度 指标 的 下 降 方向 。 其 次 ,采用 Wolfe- 
Powell 准则 自 适应 改变 控制 因子 的 大 小 ,克服 了 计算 
过 程 中 的 数值 振荡 和 不 收敛 问题 。 最 后 ,为 了 获得 随 
机 天 数 的 上 下 界 ,基于 超 参数 凸 模型 构造 了 相应 的 显 
式 渤 代 公式 。 此 外 ,通过 一 个 数学 算 例 、 个 有 限 元 
算 例 和 一 个 工业 机 器 人 应 用 验证 了 自 适应 共生 控 制 
训 涛 在 处 理 混合 可 靠 性 问题 的 有 效 性 。 结 果 表明 ,对 
扰 高 非 线性 问题 , 自 适应 共 e 控 制 方法 可 以 有 效 地 避 
件 清 算 过 程 中 的 振荡 现象 ,满足 了 计算 精度 要 求 的 同 
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